July 1994 HU Berlin-IEP-94/11 



^ 






Compact Lattice QED with Staggered Fermions 



§^ : and Chiral Symmetry Breaking □ 



A. Hoferichter, V.K. Mitrjushkin f\, M. Miiller-Preussker 



Q^ ' Humholdt-Universitdt, Institut fiir Physik, 10099 Berlin, Germany 

o' 

O 

o 

■^ ■ Abstract 

iY' Different formulations of the Ad compact lattice QED with staggered fermions 

^ ■ (standard Wilson and modified by suppression of lattice artifacts) are investigated 

by Monte Carlo simulations within the quenched approximation. We show that 

k> i after suppressing lattice artifacts the system undergoes a phase transition from the 

Vh I Coulomb phase into a presumably weakly chirally broken phase only at (unphysical) 

negative /3-values. 



i) The lattice formulation of any field theory is not unique. On a physical ground 
one has to decide which version is realized in nature, if different lattice versions do 
not belong to the same universality class. 

At present lattice QED is mostly discussed within a non-compact realization 
of the gauge part in the lattice action |jl|, |^ . However, considering QED as arising 
from a subgroup of a non-abelian (e.g. grand unified) gauge theory we are led to a 
compact description of the local gauge symmetry. 

The aim of this note is to shed some light on the mechanism of spontaneous 
chiral symmetry breaking (SCSB) in lattice gauge theories. Compact lattice QED 
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with staggered fermions serves very well for this purpose. It is comparatively simple, 
and in the limit of vanishing fermion mass the classical action is chirally symmetric. 
In case of the standard Wilson gauge action the quantized theory undergoes a phase 
transition, such that in the zero-mass limit the chiral condensate (xx) becomes 
zero within the Coulomb phase (i.e. at /5 = l/e^^^.^ > Po — 1-0 ), and non-zero 
in the confinement phase { P < Po ) 0. 

The confinement phase of compact QED does not seem to be realized in nature. 
It is related to the existence and condensation of De Grand-Toussaint monopoles 0] 
which we can understand to be artifacts of the realization of the lattice discretiza- 
tion. 

In case of the pure gauge theory it has been shown [|], ^ that a complete 
suppression of these monopoles removes the phase transition at Pq . One is left 
with a unique Coulomb phase (at least at positive /5's) 0. A further modification 
of the theory by suppressing additionally negative plaquettes was shown not to 
change the spectrum of gauge degrees of freedom but to improve the overlap of the 
plaquette operator with the lowest state, at least, in the case of the photon quantum 
numbers [|]. 

A comparative study of the chiral properties of the standard and of the modified 
theories with compact gauge fields coupled to staggered fermions is the aim of this 
investigation. Within the quenched approximation we are going to discuss the phase 
structure of these theories. We will show that the modified theories undergo phase 
transitions at negative (i. e. unphysical) /5-values only. 



ii) We consider compact f/(l) gauge fields coupled to staggered fermions defined 
on a lattice of the size L^ ■ Lt (different time and spatial extensions are only 
introduced for the computation of correlators). The gauge fields are treated with 
the Wilson action modified by suppressing lattice artifacts, such that the gauge 
action in general reads as follows 

Scm =P-Y.{l-cosep) + \K-Y. \Kc\ + \p ■ E(l - sign(cos^p)) , (1) 

P c p 

where Ui = Ux^j. = exp{i6x^), 6x^ G (— vr, -n] are the field variables defined on the 
links / = (x,/i) . The plaquette angles are 6p = O^-^y = 6x;f, + 6x+fi.y-6x+i>;ti-dx;u ■ 
The second and the third term in the action have been added in order to suppress 
monopoles and negative plaquette values, respectively, as lattice artifacts. The 
monopole currents Kc are defined as follows. The plaquette angle can be splitted 

ep = ep + 27mp, -71 < Op < TT, np = 0, ±1, ±2, (2) 

where 6p describes the (gauge- invariant) 'electromagnetic' flux through the pla- 
quette and np is the number of 'Dirac' strings passing through it. The net number 
of Dirac strings passing through the surface of an elementary 3d cube c^^^ and 
taken with the correct mutual orientation determines the monopole charge within 



this cube, i.e., a monopole current along the corresponding dual link labeled by 

(x,/i) 

Kc = K.,,^ = ttE^p = T^E(^P-2^M = -J2^P- (3) 

Ztc p 2n p p 

In the action (|l]) the parameters A^ and Xp play the role of chemical potentials 
for monopoles and negative plaquettes, respectively. In the present investigations 
we have chosen the following cases (A^, Ap) = (0, 0) - standard Wilson action 
WA, (oo, 0) - modified action MAI with total monopole suppression, and (oo, oo) 
- modified action MA2 with simultaneous total monopole and negative plaquette 
suppression. In practice we realize this total suppression via Kronecker (5's in the 
measure of the functional integral. They are easily locally taken into account in the 
heat bath or Metropolis Monte Carlo updating procedure. 

The staggered fermion part of the action looks as follows. 

Sf{Ui,x,x) = -xix)iM+m)^yXiy) (4) 

M:,y = --Y.^^(^)P^A,^+f^-Ul^^y,x-ii] (5) 

where 77^(0;) = (^^iyi+x2+-+x^-i^ rii{x) = 1 , and x ? X represent one- 
component Grassmanian variables. 

Gauge invariance requires the gauge fields to couple to fermions via the phase 
factors Uxfi , i.e., in a compact way. This holds independently of the compact or 
non-compact realization of the gauge field action So ■ Therefore, one should keep 
in mind that lattice QED within the standard compact and the non-compact case 
agree at strong coupling /? = . On the contrary, the constraints due to the 
modification terms in MAI and MA2 change the strong coupling behaviour. There- 
fore, from the beginning at strong coupling we cannot a priori expect qualitative 
agreement with non-compact QED. 

Here, we treat the fermions within the quenched approximation. We expect that 
the chiral properties of the full theory will be qualitatively correctly described. 

First of all we are going to discuss bulk properties, i.e., we measure the av- 
erage plaquette action Ep = (1 — cos 9p) as well as the chiral condensate 
ixx) = "7 ( Tr(7Vl + m)^ ) . The latter is computed applying the conjugate 
gradient algorithm combined with the noisy estimator method. 

Moreover, in order to distinguish the Coulomb phase for the modified actions 
even in the negative /3-range we calculated the J"^*" = l"*"" plaquette-plaquette 
correlator, which for non-vanishing momenta Pi = 27cki/Ls with ki =±1,±2,... 
contains the photon-state ( 1 ) as the leading contribution. The photon correla- 
tor r(r) is 

T{r) = {<!>;{t + r).%it)),, (6) 



where the operators $p'(t) are defined as Fourier transforms of sin 6p P, ^ . This 
correlator is expected to behave as 



rxr) 
r(o) 



A- 



exp(-T-E{p))+exp(-{Lt-T)-E{p)) +..., < A < 1 (7) 



with the energy E{p) in units of the lattice spacing given by the lattice dispersion 
relation 

sinh^ -^ = J2 sin^ ^ (8) 

for a zero-mass excitation. 



iii) First of all we have studied the phase structure of the two modified theories 
(MAI and MA2) in comparison with the standard theory (WA). The simplest way 
is to consider the average plaquette Ep as a function of P . We extend the range 
to be investigated to negative (i.e., a priori unphysical) P values. 

For the standard theory the plaquette action shows two transition points at 
P ~ ± 1.0 , related to the symmetry /3 ^ — (3 and Ux^ — > V^l{^) ' Ux^j. 
(full circles in Fig. la). In the same figure we have plotted (open squares) the result 
for the modified action MAI. Obviously the transition at positive (3 is removed 
leaving us only with a rather strong discontinuity at (3 ~ — 0.7 . For the 
further modification by suppressing additionally negative plaquettes (MA2) this 
discontinuity is shifted even further down to a /?c — — 1-68 (full squares in 
Fig. la). There we have seen a clear bi-stable behaviour indicating the existence of 
a first order phase transition (Fig. lb). 

The interesting lesson from these calculations is that in the strong coupling 
region (as well as at negative /?'s) there is a very strong infiuence of the lattice 
artifacts considered. 

We can ask, whether the negative /3-range /3c < /3 < can be interpreted 
as an (analytical) continuation of the Coulomb phase at positive /3 . In order to 
answer this question we computed the photon propagator as explained before. The 
results are shown in Fig. 2. For the case MA2 we have plotted r(r) for /3 = 0.9 
and P = — 1.1 , for comparison, at lowest non- vanishing momentum. The 
theoretically expected zero-mass curve is also shown (broken line). All the data 
points are in perfect agreement with a zero-mass photon. 

In Fig. 3 we have plotted the chiral condensate (xx) fo^' fixed fermion mass 
m = 0.04 (in units of the lattice spacing) as a function of P for the standard 
action WA and for the modified one MA2. Whereas for WA {xx) shows clearly 
chiral symmetry breaking at /3 ~ 1.0 (as well as a sign of the symmetric transition 
at /3 ~ — 1.0 ) it stays small and approximately constant for the modified action 
down to /3c — — 1.68 . Thus, it does not differ from the behaviour at P > 1.0 , 
i.e., in the Coulomb phase. However, we observed fiuctuations of the chiral order 



parameter becoming stronger the nearer we came to the phase transition point f3c ■ 
At (3 < Pc the fluctuations are very strong making it hard to measure (xx) ^^ 
the zero-mass hmit. The rise of the chiral condensate in this range, nevertheless, 
could indicate the possibility of a weak spontaneous chiral symmetry breaking. 

At several (3 values we have studied the mass dependence of the chiral con- 
densate in the range m = 0.01 . . . 0.05 for different lattice sizes up to 16^ (see 
Table). As an example in Fig. 4 we have plotted the mass dependence for f3 = 0. 
The finite size effects turn out to be small. The extrapolation down to m = 
provides a chiral condensate value compatible with zero. An analogous behaviour 
has been seen at /3 = 0.5, 1.1 and —1.1 . 



iv) Our main conclusion is the following. At real bare coupling, i.e., for positive 
P = 1/e^^^g the modified theory MA2 (at least in quenched approximation) is 
fully in the Coulomb phase, and there is no indication for a phase transition and 
for spontaneous chiral symmetry breaking. Therefore, we conclude that SCSB in 
the standard theory is due to lattice artifacts. This is similar to the situation of 



compact QED with Wilson fermions 0, 10 



One can try an analytical continuation to negative /? values. We have seen the 
Coulomb phase to be extended down to a certain value /3c < depending on the 
concrete realization of the modified gauge action (called MAI or MA2). There, we 
observed a phase transition presumably of first order. The phase beyond jSc seems 
to exhibit weak chiral symmetry breaking. But, more thorough investigations with 
very high statistics are required before drawing any conclusions in this respect. 

The question arises if the inclusion of dynamical fermions could change the 
picture. However, presumably one should not expect qualitative changes as we 
know from studies of both compact and non-compact QED cases. 
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Table 



The chiral condensate (xx) for MA2 at different values of /3, m and for different 
lattice sizes. The given errors are purely statistical. 



p 


L 


m = 0.01 


m = 0.02 


m = 0.03 


m = 0.04 


m = 0.05 




0.0 


8 


.0099(6) 


.0198(8) 


.0299(7) 


.0378(6) 


.04601(8) 


0.0 


12 


.00956(5) 


.0213(5) 


.0286(2) 


.0380(4) 


.0469(2) 


0.0 


16 


.00934(2) 


.01923(9) 


.0288(1) 


.03748(7) 


.04644(6) 




0.5 


8 


.00854(2) 


.01705(4) 


.02554(3) 


.03432(8) 


.04264(6) 


0.5 


12 


.0090(1) 


.01739(9) 


.0277(4) 


.0360(3) 


.0430(2) 


0.5 


16 


.00871(4) 


.0181(2) 


.0262(1) 


.0345(1) 


.0450(3) 




1.1 


8 


.00784(4) 


.0159(1) 


.0239(2) 


.0328(3) 


.0404(6) 


1.1 


12 


.008124(6) 


.0169(3) 


.0257(4) 


.0316(3) 


.0394(2) 


1.1 


16 


.0085(2) 


.0171(2) 


.0245(2) 


.0331(2) 


.0417(2) 




-1.1 


12 


.01260(4) 


.02376(3) 


.03595(4) 


.04651(4) 


.05967(6) 


-1.1 


16 


.01202(1) 


.02384(2) 


.03556(4) 


.04729(4) 


.05879(4) 




-2.0 


12 


.0343(2) 


- 


.0852(3) 


- 


.1252(2) 


-2.0 


16 


.02822(4) 


.0585(1) 


.0807(1) 


- 


.1236(1) 



Figure Captions 

Fig. la: Plaquette energy Ep as a function of f3 for standard Wilson action 

WA and modified actions MAI, MA2. Lattice sizes are 6^ for WA and MA2 and 

4^ for MAI. 

Fig. lb: Time history for the plaquette energy Ep at the transition 

point (3c = — 1.68 for the modified action MA2. The lattice size is 6^ . 

Fig. 2: The photon correlator T{t) measured at (3 = 0.9 and —1.1 with lowest 

non- vanishing momentum on a 12 ■ 6^ - lattice for MA2. 

Fig. 3: Chiral order parameter (xx) ^^ ^ function of (3 for fixed fermion mass 

m = 0.04 for WA and MA2, respectively. The lattice size is 6^ . 

Fig. 4: Mass dependence of (xx) fo^^ the modified action MA2 at /5 = 0. and 

different lattice sizes ( 8"^, 12^, 16'' ). 
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